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'^ The Spitzer Space Telescope has undertaken the deepest ever observations of 

the 24/im sky in the ELAIS-Nl field with the MIPS instrument as part of GOODS 
^ . Science Verification observations. We present completeness corrected extragalac- 

^ ! tic source counts down to 24/im fiux densities of 20/iJy (30% completeness), a 

fr^ , factor of 10000 more sensitive than IRAS. The shape of the counts confirms 

p^ ! model predictions for a strong evolution of the infrared luminosity function be- 

tJ" ! tween redshift of and 1, and suggests a fiattening in the evolutionary parameters 

^"-"^ ' at higher redshifts. Models that fit the counts indicate that luminous infrared 

galaxies (10^^ <L(8-1000/im)< 10^^ Lq) in the redshift range 0.5< z < 2.5 con- 






o 



stitute >60% of 24/im sources seen in the fiux range 20<S,y < 1000 /iJy. At least 
85%, and possibly all of the 24/im sources have a counterpart in the IRAC 3.6/im 
c^ ! and 4.5/im near- infrared channels consistent with the expected spectral energy 

distribution of infrared luminous galaxies at moderate redshift. The similarity 
k> i between the observed mid-infrared to near-infrared fiux ratios of the Spitzer de- 

tected sources and the 15/im/HK-band fiux ratios of the ISOCAM 15/im sources 
seen in the HDF-N strongly suggests that faint 24/im sources are high-redshift 
analogs of ISOCAM 15/im sources, and that they have the potential to provide 
an evolutionary connection between the well-studied 2; ~ 3 Lyman-break galaxy 
population and the dusty starburst galaxies seen at z ~ 1. 
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1. Introduction 

Deep extragalactic surveys at mid-infrared and submillimeter wavelengths have revealed 
a population of dusty galaxies with large mid- and far-infrared luminosities at intermediate 
and high redshifts (see e.g. Franceschini et al. 2003, and references therein). The number 
counts of these objects is more than an order of magnitude higher than expected if the 
local infrared luminosity function of galaxies were constant with redshift. Model fits to 
the observed counts at these wavelengths by various groups suggest a rapid increase in the 
co-moving number density of infrared luminous galaxies between redshifts of and 1 — 1.5 
followed by a fiattening or decrease out to higher redshifts (Xu et al. 2001; Chary & Elbaz 
2001; Lagache et al. 2003; King & Rowan-Robinson 2003, ; hereafter XU, CE, LDP and 
KRR respectively). The mid-infrared bandpasses i.e. ISOCAM 15/xm and Spitzer 24/im 
filters trace rest-frame 7.7/im and 12/im at ^~1 and 2;~2 respectively where the radiation is 
dominated by the line and continuum emission from polycyclic aromatic hydrocarbons and 
very small dust grains (see e.g. Dwek et al. 1997, for a review). As a result, it has been unclear 
if the evolutionary turn over at z~l is because of the /^-correction in the ISOCAM IS^um 
filter which would selectively detect objects with strong 7.7/im PAH features at that redshift 
or because of a real change in the evolutionary parameters. In this paper, we present faint 
24//m counts from S])iteer observations of the European Large Area ISO Survey-Nl (ELAIS- 
Nl) field, evaluate their contribution to the extragalactic background light and compare 
the properties of these objects with the ISOCAM detected 15/im sources with particular 
emphasis on their near-infrared/IRAC counterparts. Throughout this paper, we refer to 
luminous infrared galaxies (LIGs) as objects with 10^^ <Lir < 10^^ L© and ultraluminous 
infrared galaxies (ULIGs) as objects with Ljr > 10^^ Lq. 



2. Observations 

The Spitzer observations of ELAIS-Nl (a = 16/i09m20s; 6 = 54°57'00"; J2000) were 
undertaken under Director's Discretionary Time as part of the Great Observatories Origins 
Deep Survey (GOODS) Science Verification program to assess the effect of source confusion 
on ultradeep Spitzer surveys. The maximum integration time amounted to 4620sec per 
detector pixel of MIPS 24/im integration and 6900sec of IRAC observations corresponding to 
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background limited 5cr sensitivities of 6/iJy/pix and 0.4/iJy/pix (3.6 & 4.5/im) respectively 
(see Fazio et al. 2004; Rieke et al. 2004, for instrument description). The output of the 
Spitzer Science Center v9.5 pipeline was post-processed after making corrections for varying 
sky background, striping (jailbar) correction for every fourth column due to a bias drift 
in the readout electronics, and distortion corrections based on a grid of 2MASS stars and 
cross-channel counterparts in the IRAC and MIPS images. The images were then resampled 
(drizzled; Fruchter & Hook 2002) onto a 1.2" pixel grid for MIPS and 0.6" for IRAC and 
sources extracted. The observations covered a total area of ~165n' with an integration time 
greater than 900s. Two rectangular strips with a combined area of SOD' have the deepest 
exposures (>4200sec). Figure 1 shows a fraction of the 24|Um image with dots representing 
the IRAC counterparts of the MIPS sources. 

For the IRAC data, roughly half the field was covered by the 3.6/im and 5.8/im channels, 
while the other half was covered by the 4.5 and S.Oyum channels. Only the 3.6 and 4.5yum 
IRAC images are considered here, which together provide the best angular resolution and 
faintest flux limits currently available for studying the counterparts to the 24yum sources. 
The source density in the IRAC channels was found to be about a factor of 4 higher than in 
the 24/im image at these depths. 

Sources catalogs were generated using SExtractor (Bertin & Arnouts 1996). SExtractor 
parameters were fine tuned from simulated pure-noise images. These images were simulated 
by generating noise-only data frames where the rms for each 30s frame was the same as 
measured in the data. Each frame was then run through the post-processing pipeline, ap- 
plying the same distortion corrections and drizzle parameters in generating a final stack as 
for the real data. In general, we find that the majority of the sources in the real data are 
indistinguishable from a point source and so hereafter, we assume that all sources are point 
sources. The PSF was measured from the real data and normalized by comparing the flux 
in an 18" radius aperture of the PSF with an equivalent SIRTF/TINYTIM (Krist, personal 
communication) generated PSF with the same FWHM. Point sources were then added on 
to the pure noise image with three different model flux distributions corresponding to the 
KRR, CE and LDP models. Figure 2 illustrates the difference in source density between the 
data and the models. Allowing for a matching radius of 2.9", we flnd that our SExtractor 
parameters result in <1% of spurious sources in our simulated images where we know the 
input catalog of sources. The flux of the extracted source was measured in a circular aper- 
ture of radius 6" and corrected upwards by a factor 1.8 to account for the wings of the PSF 
performed. We found that as a result of the large source density and the 4.7" FWHM PSF, 
using larger beams results in substantial contamination of the source flux from neighboring 
sources. 
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Completeness corrections and flux calibration were measured using a thorough Monte- 
Carlo approach. In each iteration of the process, 14 artificial sources were added to the 
original data image, and SExtractor run on the resultant new image, with detection and 
flux measurement parameters identical to hose used for the original catalogs. The number 
of sources added was deliberately kept small to minimize source crowding in the images. 
Sources present in the original image were tagged in the new image depending on the change 
in their position and/or change in their flux density. Thus, if a source in the new image had 
a position within 2.9" in the artificial image and had a flux difference of less than 5% it was 
tagged as unaffected by the artificial sources. An artifical source was regarded as detected 
if one of the untagged sources matched the input position of the artificial source to within 
2.9". Detection and flux measurement parameters were identical to those described above. 
The process was repeated a 1000 times with the 14 sources randomly generated from a flat 
flux distribution between 20^Jy and 2000^Jy. A matrix Pj^ for the output flux distribution 
of the artificial sources was generated where i is the input flux and j the output flux (e.g. 
Small et al. 1995). The nature of the Pjj matrix is such that for a particular i, the sum over 
all j is less than unity. This is the completeness correction factor which was measured to 
be 50% at an input flux value of 35//Jy. The observed catalog of sources in the real image 
was then distributed among the flux bins. The Pij matrix was renormalized such that the 
sum over i for a particular j was equal to the number of detected output sources in that 
flux bin. The completeness corrected counts in each flux bin i is then the sum over j of the 
renormalized Pij matrix. This approach was tested on the three simulated images, each with 
their own Pj^ matrix, and the input and output catalog compared to ensure the accuracy of 
the approach. 



Uncertainties in the observed counts were assumed to be Poissonian i.e. v Nj where N^ 
is the number of observed counts in the Fj flux bin. These were then propagated through 
the Pij matrix to derive the uncertainty in the completeness corrected counts. Attempts 
to quantify the systematics in the counts were undertaken. These were dominated by the 
adopted size of the beam used to measure the flux and the resultant correction. A 5% 
flux calibration error was also introduced to assess the effect on the counts. The resultant 
systematic uncertainties are marked as the hatched region on the top panel of Figure 3. 

To assess the accuracy of the source extraction at the faintest flux bins, where the 
completeness fraction is ~30%, a pure noise image was generated on a 1.2" pixel grid and 
point sources added with a flux distribution that is identical to the completeness corrected 
counts at S24^m > 20/iJy but with an extrapolation to flux densities of l/xJy using the CE 
model. Artificial sources were added and extracted from this simulated image using the 
technique described earlier. The completeness corrected counts in this simulated image are 
only marginally different from those derived from the real data and hence we conclude that 
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despite the large completeness corrections, the shape of the counts at the faint end are 
reliably determined. 



3. Galaxy Counts and Extragalactic Background Light 

Table 1 and Figure 3 shows the completeness corrected 24 /im counts in ELAIS-Nl and 
the actual number of detected sources in each flux bin. Also shown in Figure 3 are the 
model 24;um counts from XU, CE, LDP and KRR. All the model counts are too high at 
bright flux densities (824^01 > 500/iJy) with the CE model being the most deviant. This 
could be attributed to cosmic variance and poor statistics on bright objects in this field. 
However, based on the Spitzer 24yum counts presented by Papovich et al. (2004) in this issue 
which offer much better statistics at these brighter flux densities, it appears that this is a 
real feature. This is most likely because the evolutionary models of the various groups were 
fit to the published ISOCAM IS^um counts from various different fields at flux densities of 
~1 mJy which had large associated uncertainties and have since been shown to be too high 
(see e.g. Gruppioni et al. 2002; Elbaz et al. 1999). At fiux densities fainter than 400 //Jy, 
our counts indicate that both the LDP model and KRR model are inconsistent with the 
observed counts. The LDP model is about 30% too low and has a faint end slope of S~^'^ 
while the KRR model is comparably low and has a faint end slope of S~^'^^. For comparison, 
the observed counts in this field at fiux densities in the range 20— 200/iJy can be fit by a 
polynomial of the form dN/dS (arcmin'^ fiJy-^)= 92x(S yuJy)-i-^±°-i°. The XU and CE 
models are within the uncertainties in the counts in this fiux range and have S^^'^ slopes as 
well. 

To quantify the quality of the model source counts fits to the data at the faint end, we 
performed a minimum sum of absolute errors, a minimum chi-square as well as a maximum 
likelihood analysis of the observed completeness corrected counts since the uncertainty is 
dominated by Poissonian errors. The likelihood, L is calculated as: 



In L = ^ m, ln(Q) - ^ 



1=0 i=0 

The third term in the calculation of the likelihood, — Sln(?77,j!) is irrelevant and is omitted 
in Equation 1. m, is the model counts in each of n flux bins and q is the observed counts 
corrected for completeness in those flux bins. We flnd that for both the minimum absolute 
errors and chi-square technique, the CE model which includes both density and luminosity 
evolution provides the best flt to the counts at flux values fainter 500 yuJy. However, the 
maximum likelihood technique indicates that the XU model is marginally a better flt. Thus, 
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we conclude that the CE and XU models both provide reasonably good fits to the observed 
data at the faint end. Although neither of these models provide good fits at 824^111 > 500/i Jy, 
the counts at these fiux limits are dominated by galaxies at z<0.5. Roughly speaking, the 
implication of this poor fit is that at these redshifts, the LIG component is evolving much 
more rapidly than the ULIG component. However, the paucity of sources at faint 24/im 
fiuxes implies that the rapid increase in the number density of infrared luminous galaxies 
between redshifts of and 1 — 1.5 that is required by the models does not extend to higher 
redshifts and the evolution either reaches a plateau or declines. Since the optical/ultraviolet 
luminosity density (with no dust corrections) also rises but at a slower rate between z~0 
and z^l (see Steidel et al. 1999, and references therein), the implication is that extinction- 
corrections to UV-derived star-formation rate measurements are increasing with redshift, 
particularly between 0< 2; < 1. 

The predicted redshift distribution of the sources (Figure 4) in this fiux range as derived 
from the models shows a bi-modal structure, peaking at redshift 1 and redshift 2 with ~20% 
of sources being at redshifts between 2<z<3. The bi-modal nature of the redshift distribution 
is a result of the two strongest set of PAH features, those around S/xm and 12/im, being 
redshifted into the 24|Um passband. Luminous infrared galaxies account for about 60% of all 
sources while ULIGs are about 10%. The remainder consists of M82-analogs, i.e. low-redshift 
objects with L(8-1000/im)< 10^^ Lq. 

The robustness of this redshift distribution can be questioned since it depends on empir- 
ical evolutionary models which broadly fit other multi-wavelength counts but provide poor 
fits at bright 24;um fiux densities. CE emphasized the degeneracy in evolutionary parameters 
for the infrared luminosity function i.e. a number of models which provide reasonable fits to 
the observed counts can yield varying counts/star-formation rate estimates. However, the 
key conclusion was that all models required a fiattening in their (1+z)" evolution at z~l. 
This is confirmed by the Spitzer observations. We have re-derived the evolutionary parame- 
ters to fit the 24/im counts in this paper and Papovich et al. (2004) as well as the far-infrared 
counts in Dole et al. (2004). We find that the key variable is the ratio of LIGs to ULIGs 
which we now evolve independently unlike in CE. Using the new evolutionary parameters 
which we will discuss in a later paper, we find that although the fractional contribution of 
LIGs and ULIGs can differ slightly from that presented in Figure 4, the shape of the overall 
redshift distribution remains the same. 

The CE models derive the total 24yum EBL to be 3.7nW m^^ sr^^of which 32% comes 
from sources brighter than lOOOyuJy. Since the CE model is a factor of ~4 too high at the 
bright fiux densities, the 24/xm EBL must be significantly lower than this estimate. The 
24//m IGL from the detected sources in the fiux density range 20<Siy < 1000 yuJy in these 
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Spitzer ohservations is 1.8±0.2 nW m~^ sr^^ while the completeness corrected counts yield 
a value of 2.0±0.2 nW m~^ sr~^. Thus, we conclude that these deep 5]9^teer observations 
have resolved out at least 50% of the 24|Uni EBL and probably as much as 65%. 



4. Near-Infrared Counterparts of 24fim Sources 

ISOCAM 15//m sources appear to be disparate from the SCUBA 850|Um sources in 
that they have relatively bright optical/near-infrared counterparts which span the range of 
Hubble-types in morphology. This seems to suggest that they are typically massive galaxies 
with stellar mass > 10^'^ Mq (Franceschini et al. 2003). In contrast, the SCUBA 850/xm 
sources are optically faint and mostly irregular systems (Chapman et al. 2003). This could 
be a result of morphological /i;— correction since visible light observations of SCUBA sources 
trace rest-frame ultraviolet emission which could be patchy due to dust extinction internal to 
the galaxy. Since the models suggest that the 24^m sources presented here are at redshifts 
intermediate to the ISOCAM and 850/um sources, this would imply that they have rest-frame 
properties that are intermediate to those sources. 

We matched sources between MIPS and IRAC using a nearest neighbor match between 
the 24|Um catalog and the 3.6 and 4.5 yum catalogs (Figure 5). At a 24yum flux limit of 50//Jy, 
more than 85% of the sources have IRAC counterparts within 3", as compared to the 25% 
expected for random matches to objects with the same density; thus at least 60% of the 
matches are physical. 

Direct inspection of the remaining 15% shows that in almost all cases, separations of 3" 
or larger arise not from a lack of physical counterpart, but due to a variety of other factors: 
7.8% of MIPS sources fall in a region of low exposure time in either the IRAC image or MIPS 
image or both; 3.6% are cases where the IRAC counterpart was not automatically detected 
because it falls near a much brighter IRAC sources and was not deblended from the brighter 
IRAC sources; 6.5% of MIPS sources appear to be blended emission from two or more 
nearby IRAC sources in which case the MIPS source appears slightly extended/distorted 
and its centroid falls between the IRAC positions which are more than 3" away. There are 
however 0.9% of cases where the MIPS source genuinely appears to have no IRAC counterpart 
to the sensitivity of the IRAC image. Since the distribution of MIPS/IRAC flux ratios is 
broad, depending on the nature of the source and its redshift, faint MIPS sources could 
have counterparts which are below the flux threshold of the IRAC observations. Therefore 
we conclude that the fraction of 24/im sources with reliable near-infrared associations is 
-100%. 



Having identified the counterparts of tlie sources, we tlien proceed to estimate a flux 
ratio between tlie mid-infrared and near-infrared wavelengtlis. For tlie MIPS catalog we 
adopted a beam size corrected flux value as described in Section 2. For the near-infrared 
we adopted SExtractor flux values in a 6" diameter beam and applied a 20% beamsize 
correction based on estimates published by the Spitzer Science Center. Figure 5 illustrates 
the mid-infrared to near-infrared flux ratios of the Spitzer 24/im sources detected here. Also 
shown as solid circles is the 15//m to HK-band flux ratio (partly from Barger et al. 2002, 
and partly from the GOODS data) for the ISOCAM 15^m sources in the HDF-N which 
have spectroscopically determined redshifts peaking at ~0.8. Due to the shape of the dust 
spectral energy distribution at these wavelengths, sources at the same redshift will mostly 
be brighter at observed 24/im than at observed 15/im. However, at 2; ~ 2, the properties of 
24/im sources should trace the same rest-frame emission as the 15/im observations at z~l. 
Since, the mid-infrared flux of a galaxy traces its total star-formation rate, while the near- 
infrared flux is a measure of the stellar mass in the galaxies, the similarity in this flux ratio 
between the ISOCAM and MIPS sources broadly suggests that the specific star-formation 
rate i.e. the ratio of SFR to stellar mass, is very similar among the z^l ISOCAM galaxies 
and z^l-2 MIPS galaxies. 
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Fig. 1. — Image showing the quality of the MIPS 24^m data in the ELAIS-Nl field. The 
stretch is in flux values in units of microJy. The dots indicate the position of the IRAC 
counterpart. 
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Fig. 2. — Comparison between the source densities in the real data (center) and the simu- 
lated image from the King & Rowan-Robinson model (top) and the Chary & Elbaz model 
(bottom). Clearly the KRR model predicts too many faint sources while the CE model has 
too many bright sources with respect to the data. 
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Fig. 3. — Completeness corrected galaxy counts in the MIPS 24/ini channel from Spitzer 
observations of ELAIS-Nl. The error bars reflect the Poissonian uncertainty. The horizontal 
bars represent the minimum and maximum flux density in that bin. The lines show four 
preferred models for 24/im counts (see text). The symbols are plotted at the average of the 
flux densities of the detected sources in that bin for the data while the hues are plotted 
at the counts-weighted flux average for the models. The lower plot in the figure shows the 
histogram of the actual number of sources detected in each flux bin without any completeness 
correction. 
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Fig. 4. — Redshift distribution expected for 24/im sources in the flux range 20< S",^ < 500/iJy 
as derived from the best fitting CE modeL The dotted line indicates the contribution from 
ULIGs (LiR =L(8-1000;um)> 10^^ L©), the dashed hne is for LIGs with the triple-dot-dashed 
hne the contribution from low-luminosity starbursts with Lir < 10^^ Lq. The size of the 
redshift bins is 0.1. 
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Fig. 5. — Ratio between the mid-infrared (MIR; 24^m) flux density and the near-infrared 
(NIR) flux density in the 3.6^ni (plus) or 4.5/im (triangles) IRAC channels for all 24/im 
sources brighter than SOyuJy. Also shown in the lower plot is the distribution of separation 
in arcseconds between the MIPS source and its IRAC counterpart. Of the ~1000 sources 
brighter than 50/xJy in our 24yum catalog, more than 85% have associations in the IRAC 
images which are less than 3" away. The plot also shows the IS^um to HK-band flux ratio for 
the ISOCAM sources in HDF-N (solid circles). To qualify the nature of the MIPS sources, 
we also show the evolution of the colors of two nearby galaxies; the ultraluminous starburst 
Arp220 (solid squares) and the quiescent starburst/spiral galaxy M51 (solid triangle) as a 
function of redshift. The redshift increases in steps of 0.25 with the redshift of the first point 
as labelled. Arp220 and M51 at lower redshift will be to the right of this plot. 
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Table 1. 24/im Counts in the ELAIS-Nl field'^ 



^low 


"^high 


AvgS. 


Observed dN/dS 


Corrected dN/dS 


Uncertainty dN/dS 


20.0 


30.0 


24.8 


1.2E-01 


3.6E-01 


2.5E-02 


30.0 


45.0 


37.1 


8.9E-02 


2.0E-01 


1.3E-02 


45.0 


67.5 


55.2 


5.9E-02 


l.OE-01 


6.8E-03 


67.5 


101.2 


83.8 


3.8E-02 


5.3E-02 


3.6E-03 


101.2 


151.9 


123.1 


2.4E-02 


2.8E-02 


2.0E-03 


151.9 


227.8 


185.9 


1.5E-02 


1.5E-02 


l.lE-03 


227.8 


341.7 


277.9 


6.8E-03 


6.4E-03 


5.7E-04 


341.7 


512.6 


415.8 


3.0E-03 


2.5E-03 


2.7E-04 


512.6 


768.9 


597.9 


5.6E-04 


4.9E-04 


l.OE-04 


768.9 


1153.3 


891.2 


1.3E-04 


1.2E-04 


4.3E-05 



^Su is in /ijy while dN/dS is in # arcrnin ^ /ijy ^. Uncertainty in dN/dS is the Poissonian 
noise propagated through the P^j matrix. 



